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ABSTRACT 

X-ray timing of the accretion-powered pulsations during the 2003 outburst of the accreting miUisec- 
ond pulsar XTE J1814-338 has revealed variation in the pulse time of arrival residuals. These can 
be interpreted in several ways, including spin-down and wandering of the fuel impact point around 
the magnetic pole. In this Letter we show that the burst oscillations of this source are coherent with 
the persistent pulsations, to the level where they track all of the observed fluctuations. Only one 
burst, which occurs at the lowest accretion rates, shows a significant phase offset. We discuss what 
might lead to such rigid phase-locking between the modulations in the accretion and thermonuclear 
burst emission, and consider the implications for spin variation and the burst oscillation mechanism. 
Wandering of the fuel impact hot spot around a fixed magnetic pole seems the most likely cause for 
the accretion-powered pulse phase variations. This means that the burst asymmetry is coupled to the 
hot spot, not the magnetic pole. If premature ignition at this point (due to higher local temperatures) 
triggers a burning front that stalls before spreading over the entire surface, the resulting localized 
nuclear hot spot may explain the unusual burst and burst oscillation properties of this source. 
Subject headings: binaries: general, stars: individual (XTE J1814-338), stars: neutron, stars: rotation. 
X-rays: bursts. X-rays: stars 



1. INTRODUCTION 

The accreting millisecond X-ray pulsars (AMXPs) are 
a small class of neutron stars in Low Mass X-ray Bi- 
naries that show pulsations in outburst, thought to be 
caused by magnetic channeling of accreting plasma. De- 
tailed timing studies of these stars reveal diverse be- 
havior that can be interpreted in terms of spin varia- 
tion or shifts in emission pattern. Both processes may 
play a role, and the degree to which we can be confi- 
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Some additional way of verifying the timing analysis 
obtained from the accretion-powered pulsations would 
be useful, and in this respect two of the AMXPs are par- 
ticularly valuable. SAX J1808.4-3658 (J1808) and XTE 
J1814-338 (J1814) also show thermonuclear-powered pul- 
sations, or burst oscillations. These are high frequency 
variations seen during Type I X-ray bursts, powered by 
unstable burning of accreted fuel. In these systems the 
burst oscillation frequency is at (J1814), or very close 
to (J1808), the spin fre quency (|Chakrabartv et a"l]|2003l : 
IStrohmaver et aDl2003l . hereafter S03). The frequency is 
stable in the decaying tails of the bursts, with no sign of 
the larg e frequency drifts s een in other burst oscillation 
sources (|Muno et al.ll200^ . In J1814 the frequency is 
also stable during the rising phase of the bursts, making 
it the most straightforward ca ndidate for burst oscilla- 
tion timing (jWatts et al.ll2005l . hereafter W05). 

J1814's accretion-powered pulsations show significant 
pulse time of arrival (TOA) r esiduals even after c orrec- 
tion for orbital Doppler shifts (jPapitto et al.l l200^ here- 
after P07). The cause is still a matter of debate: P07 

^ A.L.Watts@uva.nI 



interpreted the observations in terms of a steady spin- 
down coupled with some jitter due to wandering of the 
fuel impact hotspot around the magnetic pole. However 
there are other possibilities, such as changes in beam- 
ing due to the accretion shock, that may also lead to 
the observed variation. In this respect analysis of the 
burst oscillations may be simpler: although the process 
responsible is not yet understood, the thermal spectrum 
suggests a purely surface mechanism, with the accretion 
shock contributing little to the observed asymmetry. 

Some initial investigation of this issue was carried out 
by S03, who found that burst oscillations in the first 12 
bursts were phase-locked to within 2.5° of the persis- 
tent pulsations. Their analysis, however, covered only 
the first 10 days of the w 50 day outburst; before any of 
the variation reported by P07 is apparent. The level of 
coherence between the two types of pulsations is also im- 
portant in our efforts to understand the burst oscillation 
mechanis m, which remains mysterious (see the r eviews 
by Stroh maver fc BildstenI l2006l iGallowav etldl [20061. 
The AMXPs are the only sources in which we can quan- 
tify the role of the magnetic field. We want to know, 
for example, whether the burst oscillations couple to the 
magnetic field or to the fuel stream impact point. 

2. TIMING ANALYSIS 

J1814 was discovered in 2003 in the Rossi X-ray Tim- 
ing Explorer (RXTE) Galactic bulge monitoring cam- 
paign (jMarkwardt fc Swankll2003( ). and remained in out- 
burst for nearly 2 months. The pulsar has a spin fre- 
quency of 314.4 Hz and resides in a bina ry with an or- 
bit of 4.3 hours (jMarkwardt et al.ll2003[ ). During the 
outburst over 425 ks of high time resolution data were 
taken with RXTE' s Proportional Counter Array (PCA, 
iJahoda et al.l[2006[ ). A total of 28 X-ray bursts were de- 
tected, all with burst oscillations at the spin frequency. 



2 



Watts, Patruno & van der Klis 



Both accretion-powered pulsations and burst oscillations 
have a strong overtone at twice the spin frequency (803). 

For our timing analysis we use all available pointed ob- 
servations from the RXTE PCA with Event mode data 
(time resolution 122 /is, 64 binned energy channels) or 
Good Xenon mode data (time resolution 1 /^s, 256 un- 
binned energy channels). The latter were rebinned in 
time to 122 /is time resolution. Data were barycen- 
tered using the JPL DE405 ephemeris and a spacecraft 
ephemeris including fine clock corrections whi ch together 
provi de an absolute timing accurac y of 3-4 ^s fiRots et al.l 
using the source position of lKrauss et al.l ()2005[ ). 

In analysing the accretion-powered pulsations we dis- 
card the X-ray bursts, removing all data from 50 s before 
to 200 s after the burst rise, and select only photons in 
the 2.5-17 keV range to maximize the signal to noise ra- 
tio. Pulse profiles are built using the fixed frequency 
solution of POT to fold segments of approximately 500 s 
of data after subtracting the background contamination 
(using the FTOOL pcabackest). The TOAs were then ob- 
tained by cross-correlating the folded profiles with a pure 
sinusoid whose frequency represents the spin frequency 
of the neutron star. The same procedure is repeated for 
the first overtone. The fiducial point used in measuring 
the TOAs was the peak of the sine wave being cross- 
correlated. The determination of pulse TOAs and their 
statistical uncertainties follows the standard radio pulsar 
technique (|Tavloij|1993D . Fitting a Keplerian orbit plus a 
constant spin frequency or spin derivative v, we obtain 
solutions consistent with those of P07. 

Like P07, we find a strong anti-correlation between 
flux and the accretion-powered pulse TOA residuals. 
It is interesting, however, that this anti-correlation be- 
comes weaker when one considers the residuals from an 
ephemeris that includes a constant ly term. A rank corre- 
lation test between flux and residuals for a constant fre- 
quency model gives a Spearman coefficient oi p ~ —0.71 
(fundamental) and p — —0.90 (first overtone), with a 
probability of < 0.01% that the two variables are not 
anti-correlated. Including spin-down, the magnitudes of 
the Spearman cocfhcicnts fall to p — —0.56 (probabil- 
ity still < 0.01% that the two variables are not anti- 
correlated, but larger than the probabilities for the zero 
spin derivative case) and p = —0.12 (a probability of 27% 
that the two variables are not anti-correlated) respec- 
tively. As we will argue in Section [3l an accretion-rate 
dependent hot spot location may be able to explain the 
entire residual record, with no need for spin variation. 

We then apply the san ie timing procedure to the X -ray 
bursts. As in W05 and IWatts fc Strohmaval (|2006f ) we 
use data where the count rate is at least twice the pre- 
burst level. For the first 27 bursts there is no evidence 
of frequency variability during the bursts (W05), so we 
use the entire burst to generate a folded profile. The fi- 
nal burst requires more care, as there is a statistically 
significant frequency drop in the late stages of the burst 
rise (Fig. 20 of W05). For this burst we use only the first 
2 s of the burst rise, before the frequency starts to shift. 
For each folded burst profile we then compute residuals 
using the same ephemeris that we used for the accretion- 
powered pulsations. Again we cross-correlate the folded 
profile using both a fundamental and an overtone. Fig- 
ure [1] shows the TOA residuals for the accretion-powered 
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Fig. 1. — Phase residuals for the accretion-powered pulsations 
(bars; red) and the burst oscillations (crosses; blue), compared to 
a model that has a constant frequency and orbital Doppler shifts. 
The final burst, which has no detectable 1st overtone component, 
requires special treatment; see text. The units on the right axis are 
rotational cycles. The top panel shows the outburst light curve. 



pulsations and the burst oscillations. 

To check whether the TOAs of the burst oscillations are 
consistent with having the same temporal dependence as 
the accretion-powered pulsations we fit a constant fre- 
quency model to the two TOA sets separately (exclud- 
ing the final burst, see below). The fitted frequency is 
the same within the statistical uncertainties both for the 
fundamental and the first overtone. We also tried fitting 
the two TOA sets with a spin frequency plus a frequency 
derivative. Again the two solutions are consistent within 
the statistical uncertainties. However P07 already noted 
that both ephemerides are a poor fit to the accretion- 
powered pulse TOA residuals, and the same is true for 
the burst oscillation residuals. They are formally incon- 
sistent with a constant frequency model: for the funda- 
mental this assumption gives a of 459 for 23 degrees 
of freedom (dof). They are also inconsistent with a con- 
stant i>, as ill this case is still large, 80 for 22 dof. 
Similar results arc obtained for the first overtone. 

It is clear from the Figure that, for all but the fi- 
nal burst, the burst oscillation TOA residuals track the 
accretion pulsation residuals. To test the phase- lock 
we computed the phase difference A0m — t/[TOAbur — 
TOAacc], using the accretion pulse TOA residuals from 
the observation containing each burst. We then fitted a 
constant. For the fundamental we found a of 15.3 
(26 dof), with best fit = (0.004 ± 0.002) rota- 

tional cycles. For the first overtone we found a of 
11.3 (26 dof), with best fit Aipm = -(0.001 ± 0.003) ro- 
tational cycles. The fact that a constant Aipm is such 
a good fit confirms that the two sets of pulsations are 
phase-locked. The small {2a) non-zero offset in the fun- 
damental bears comment. Unlike the burst oscillations, 
the accretion-powered pulsations have soft lags across 
the 2.5-20 keV band, with high er energy photon s arriv- 
ing carhcr in phase (| Watts fc Strohmaveill2006l ). Tak- 
ing this into account, we find that the burst oscillation 
TOAs are completely coincident with the softer (2.5 - 
5 keV) component of the accretion pulsations, thought 
to originate from stellar surface (Section [3]). The only 
exception to the phase-locking is the final burst, with 
A(pm = (0.2 ± 0.04) rotational cycles. 

When calculating burst oscillation residuals some care 
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is required, since if accretion continues during the burst, 
there might stiU be a contribution from the accretion- 
powered pulsations. The resulting bias in the burst os- 
cillation phase can be calculated easily by considering 
the profile that results from the addition of two oH'- 
set sinusoidal profiles. Standard trigonometric identi- 
ties yield a relation between Ayj^ (the measured offset 
between burst oscillation phase and accretion-powered 
pulse phase) and A(^;, (the bias, i.e. the offset between 
measured and true burst oscillation phase caused by 
residual accretion): 

sin(A.^™ - A(pb) 

tanAl^m — 7T7 7T7 ] 7T 7 7- 

[A'acc?'acc/A'bur?'burJ + C0S(A(y5„i - Al/Jfc ) 

where iVacc and A^bm- are the number of accretion and 
burst photons in the folded profile, race and rbur being 
the associated fractional amplitudes. The quantity A(p„j 
was measured earlier: Aipm ^ 0.01 cycles at the 99% 
confidence level. Using the values of iVacc, Nhui, '"acc 
and rbur from Table 1 of W05, Equation ^ gives a 99% 
confidence upper limit on the bias introduced by any 
residual accretion pulsation of 10~^ cycles. This is suffi- 
ciently small that it does not affect our analysis. Similar 
conclusions can be reached for the overtone. Note that 
in computing these limits we assume that the accretion 
flow parameters (A'^acc, face, A(/9m) are unchanged during 
a burst. If accretion is inhibited during a burst, due to 
radiation pressure, the bias will be lower. 

3. DISCUSSION 

J1814 has unusual burst oscillation properties com- 
pared to other sources: they occur in hydrogen- 
rich bursts, have negligible frequency and amplitude 
variation, and have a soft spectrum (S03, W05, 
IWatts fc Strohmaveil l2006f ). Our analysis has shown 
that the burst oscillations are also phase-locked to the 
accretion-powered pulsations (to within 3° at 99% con- 
fidence) despite the substantial phase-wander exhibited 
by the latter over the course of the outburst. 

In fact the burst oscillations are not only phase-locked 
but also coincident, having the same phase as the soft 
(2.5-5 keV), lagging, part of the accretion-powered pul- 
sations (although they are also at the 2a level coincident 
with the entire 2.5-20 keV accretion-powered pulse). De- 
tailed modeling of the accretion pulsations has yet to be 
done for J1814, but modeling for other AMXPs with sim- 
ilar pulse properties suggests that the soft pulsed com- 
ponent comes from a hot spot on the stellar surface, with 
the hard com ponent originating in the accretion fun- 
nel (|Poutanen fc G ierlihski 2003: iGierlihski fc PoutanenI 
120051: iFalanga fc Titarchuk,2007l i. 

We first consider what the exceptional degree of phase- 
locking implies about the cause of the variation in the 
TOA residuals. There are several parts of the system 
whose variation might affect both types of pulsation: 
surface rotation, the accretion funnel/disk, the magnetic 
field, and the fuel deposition footprint. 

Case 1: Genuine changes in the spin rate of the stellar 
surface. Our result would be consistent with a model 
where all of the variation is due to spin changes, both 
sets of pulsations being locked to the surface. How- 
ever this requires alternating spin-up and spin-down with 



\v\ ^ lO^^'^ Hz/s. Even if the crust were decoupled from 
the fluid core, this is high compared to what is achiev- 
able fro m the expected accretion or gravitational wave 
torques (jAndersson et a l. 2005; Bildstcn 1998a). Fitting 
a constant spin-down term v k Qx 10^^^ Hz/s, as argued 
by P07, does improve the quality of the fits somewhat. 
However spin derivatives ~ 10"^^ Hz/s would still 
be required on shorter timescales to explain the remain- 
ing excursions. Our results would also be consis tent with 
precession, but modeling bv Chung et al.l ()2008l ) suggests 
that precession is unlikely in this source. 

Case 2: Changes in beaming /scattering by the accre- 
tion funnel or disk. The accretion shock in the funnel 
is thought to contribute to the pulsed emission of the 
accretion-powered pulsations, leaving a signature of hard 
emission in the spectrum. If the funnel were to have 
a similar effect on the much stronger burst emission it 
would hav e to do this without lea v ing any trace in th e 
spectrum (jStrohmaver et al] 120031 : iKrauss et al.ll2005h . 
This does not seem feasible. Our result also casts doubt 
on the accretion disk being the source of the soft lag- 
ging component of the accretion-po wered pulsations (one 
of th e possibilities considered by iFalanga fc Titarchuld 
l2007f ). since it is hard to understand why the burst os- 
cillations (a surface process) would track a component 
generated in the surrounding environs. 

Case 3: Wander of the magnetic pole, or changes in 
field geometry. Motion of the magnetic pole would af- 
fect location of the accretion hot spot. If the magnetic 
field also determines the location of the nuclear burn- 
ing hot spot through modulation of ignition or emis- 
sion, this could also explain our result. However, the 
observed variability would require localized burial or am- 
plification of the poloidal field component on timescales 
of order a day. The accretion rate in the peak of 
the outb urst is at most a few percent of the Edding- 
ton rate (jGallowav et al.l [20041 ). Current modehng sug- 
gests that this is insufficient to cause burial of the polar 
field on the required t i mescales ([Brown fc Bildsten|[l998[ : 
ICumming eral1l2001l : iPavne fc Melatosll2004f ). There is 
also no obvious mechanism for field amplification: mate- 
rial arriving via a funnel flow will have almost no angu- 
lar momentum differential compared to the stellar sur- 
face. The heating associat ed with accretioii could bring a 
buried field to the surface ('Cumming et al.ll2001f ). but as 
previously stated burial is unlikely at such low accretion 
rates. 

Case 4- Wander of the fuel deposition point around 
the magnetic pole. Simulations of funneled accretion 
have shown that the fuel deposition point can exhibit 
phase excursions around a fixed magnetic pole as ac- 
cretion rate fluctuates, particularly for small misalign- 
ment angles between the magnetic pole and the spin axis 
([Romanova et al.ll2003L [20041 : iLamb et al.ll2008^ ^ Such a 
model might neatly account for the correlation between 

^ The pulse profile modeling for this source which has 
been attempted has suggested large mi salignment angles 
lIBhattacharvva et al.l l2005l : ILeahv et aU[20M) . If this were the 
case, the observed phase variability would require the fuel impact 
point to migrate back and forth by several km over the course of 
the outburst - an uncomfortably large amount. It seems more likely 
that some of the assumptions in these models, particularly the use 
of single temperature circular hotspots (known to be problematic, 
seC iWatts fc Strohmaveni2006i ). need to be revisited. 
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residuals and flux without requiring any non-zero since 
the stable position for the fuel impact point will vary 
in azimuth depending on accretion rate. The fact that 
the anti-correlation between flux and residuals is stronger 
when we set v — ^ supports this idea. 

If the most plausible explanation for the variability in 
the two sets of pulsations is the last, what physical mech- 
anisms might lead to phase-locking between the fuel im- 
pact point - which moves with accretion rate relative to 
field geometry on timescales of a few days - and the nu- 
clear burning hot spots? 

One possibility is some degree of magnetic confine- 
ment, leading to accumulation of fuel at the accretion 
hot spot. Material deposited near the polar cap will be 
prevented from spreading until the over-p ressure is suffi- 
cient to distort the field lines (jBrown fc Bi ldsten 1998,), 
even if the impact point is not precisely on the polar cap. 
To ensure that the magnetic propeller effect does not dis- 
ru pt accretion, the magnetic field of J1814 must be < 10 ^ 
G (|Psaltis fc Chakrabartvlll999tlRaT)paport et al.ll2004D . 
This is consistent with the upper limit on spin-down in- 
ferred from pulse timing (P07). At this upper limit he- 
lium could be confined until a column depth ~ 10^ g 
cm~^, but this is still well before the material re aches the 
bursting layer at column depth ~ 10^ g cm~^ (Bildsten| 
[T99 8b). Hydrogen will spread even more easily. Based 
on these estimates we conclude that fuel confinement is 
not effective, although it has been advanced as a pos- 
sible explanation for the short burst recurrence times 
(iCxallowav et al.l[2?)0l . 

The other factor distinguishing the fuel impact point 
is its temperature, which is higher than the rest of the 
star. The magnitude of the temperature differential has 
not been determined observationally, but it could cer- 
tainly affect burst emission. One possibility is that the 
higher temperatures modify the local composition (via 
steady burning between bursts, for example). The higher 
starting temperature and/or composition could in princi- 
ple modify the flux from the burst once it starts; perhaps 
driving more energetic reactions or enhancing convection. 
Whether this effect would be large enough to explain the 
high fractional amplitudes is, however, not clear. 

An alternative is the effect that a higher local tem- 
perature would have on ignition conditions. Previous 



studies that concluded that ig nition would occur pre- 
dominantly near the eq uator (jSpitkovskv et al.l 120021 : 
ICooper fc Naravanl 120071 ) did not consider the effect of 
non-uniformities in temperature. A small increase in 
temperature can have a large effect on the column depth 
required for ignition, with the hotte r area requiring a 
lower column depth ()Bildstenlll998bD . Material at the 
fuel impact point could therefore reach ignition condi- 
tions well ahead of the rest of the star. In this scenario 
the burning front might stall before spreading across the 
whole star, depending on the rate of heat transfer across 
the burning front. This would result in a brightness 
asymmetry centered on the fuel impact point. Premature 
ignition at the fuel deposition point, followed by stalling, 
could explain several observational feat ures: the rather 
small black body radius of the bursts (jGalloway et al.l 
120061). the short e r tha n expected burst recurrence times 
(|Gallowav et al.l |2004[ ). and t he burst shapes, which 
suggest off-equatorial ignition (jMaurer fc Wattal [20081. 
Note that off-equatorial ignition alone is not sufficient to 
explain the presence of an asymmetry in the burst tail: 
something else, such as stalling, is required. 

Whatever the mechanism, it must fail once the accre- 
tion rate drops (or perhaps when the burst is more ener- 
getic), since the oscillations in the final burst are offset. 
The fuel impact footprint is by this time thought to be 
smaller, since the fractional amplitude of the accretion- 
powered pulsations is rising (W05). In addition the ac- 
creted fuel has more time to spread and equilibrate be- 
tween bursts. Both factors will act to reduce the temper- 
ature differential between the fuel stream impact point 
and the rest of the star, which would make ignition at 
the fuel impact point less likely. The oscillation prop- 
erties of this final burst (substantially lower fractional 
amplitude, frequency drift) are very different to the rest 
of the sample, and it is quite plausible that a different 
burst oscillation mechanism operates in this final burst. 



We thank Mariano Mendcz and the participants of the 
workshop "A Decade of Accreting Millisecond X-ray Pul- 
sars" for the lively discussions that prompted and in- 
formed this work. 
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